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Determination of constant lattice spacing trajectories in lattice QCD 

Alan C. Irving®', UKQCD Collaboration* 

^Theoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, PO Box 
147, Liverpool L69 3BX, UK. 

We argue that lattice simulations of full QCD with varying quark mass are best conducted at fixed lattice 
spacing rather than at fixed /3. We present techniques which enable this to be carried out effectively, namely the 
tuning in bare parameter space and efficient stochastic estimation of the fermion determinant. Results and tests 
of the method are presented. We discuss other applications of such techniques. 


1. CONSTANT LATTICE SPACING 

Hadron spectrum measurements by the 
UKQCD collaboration |j^,^ using an improved ac¬ 
tion have shown a surprisingly strong dependence 
of the effective lattice spacing on the bare quark 
mass. Recent preliminary measurements of the 
Sommer scale tq using the non-perturbatively im¬ 
proved value of Csw in 16^ x 32 simulations have 
confirmed this. See Table Since the parameters 


Table 1 

Lattice spacing determined from ro |§| at /3 = 5.2 
for (!7(a)-improved Wilson fermions 
~ T330 1340 l345 T350 

a [fm] .154(2) .140(4) .130(2) .1173(16) 


in chiral perturbation theory inherit this underly¬ 
ing a-dependence, this complicates chiral extrap¬ 
olations of simulation measurements and obscures 
comparisons with quenched calculations. Indeed, 
one would not attempt to compare quenched and 
full QCD at the same value of /3. Because of the 
strong and non-trivial dependence on the the bare 
quark mass, one should avoid using fixed l3 within 
the dynamical sector also. 

With this in mind, we investigate how one 
might control the effective lattice spacing by tun¬ 
ing the bare action parameters while the effects 
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of decreasing quark mass are studied. In previous 
work [^|j^ we have described how stochastic esti¬ 
mates of the fermion determinant may be used 
to tune action parameters so as to maintain fixed 
the value of some observable F. We would suggest 
that, in order to hold the effective lattice spacing 
fixed, a suitable observable is the Sommer scale 
cq. This is free from complications due to valence 
quarks (partial quenching) and chiral symmetry 
constraints. 

For small parameter changes, one may apply a 
first order cumulant expansion ^ 

< F > 2 =< F >1 -b < FAi2 >1 -b ... 
and Ai 2 = Si — S 2 , F = F— < F > etc. (1) 


Here, < .. >i means the expectation value with 
respect to action Si. For the improved Wilson 
action. 


Ai 2 = (/32 — (di)Wa + T{P2, 1 ^ 2 ) — T{(3i, ki) . (2) 

The operator T (= TrLnMfM) is /3-dependent 
only when Csw 7 ^ 0. We may use (0) to estimate 

- df 

« < F-— > and (3) 

OK 

~ - df 

« < F{Wu + —4w) > . (4) 

O'Csw 

Matching F at nearby points {/3i, ki) and {P 2 , ^ 2 ) 
thus means solving 


d<F> 

Bk 

d<F> 

dp 


5P_ 

5k 


< F6f >1 
5k < FWa >1 ’ 


5P = P 2 — Pi etc. (5) 


(in the case where Csw = 0). This in turn requires 
reliable (stochastic) estimates of T. 
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2. LANCZOS QUADRATURE 

Bai, Fahey and Golub |Q have shown how to 
estimate the trace of a matrix function via Gaus¬ 
sian quadrature using an elegant relationship be¬ 
tween orthogonal polynomials associated with the 
measure and the Lanczos recursion scheme. We 
construct a noisy estimator for TrLn(M'lM): 

N4, N 

Et = = ^w|Ln(Aj). (6) 

^ i=l j=l 

Here {A®} {j = 1,2... A) are the eigenvalues 
of the A-dimensional tridiagonal Lanczos matrix 
formed using the noise vector (j)i as a starting vec¬ 
tor. The weights {w|} are related to the corre¬ 
sponding eigenvectors 0. In fact, LOj is just the 
first component of the jth eigenvector of the tridi¬ 
agonal matrix. Rather than use bounds for the 
truncation error associated with finite values of 
A 0, we take advantage of the rapid (exponen¬ 
tial) convergence of the quadrature to obtain 6 
figure accuracy for each noise term using 

typically A = 60 to 100 Lanczos iterations. This 
is made possible by a remarkable stability to the 
round-off error and loss of orthogonality normally 
associated with a standard Lanczos analysis of 
the eigenvalue spectrum. We have noted an ‘ex¬ 
clusion principle’ associated with the quadrature 
weights uji which ensures that only genuine eigen¬ 
values of the tridiagonal system contribute to the 
quadrature even in the presence of ‘ghost’ and 
‘spurious’ eigenvalues [^. 

3. QCD MATCHING RESULTS 

Given an initial QCD simulation at {/3 i,ki) 
there are two obvious ways to locate a nearby 
point (/ 321 K 2 ) on the constant F curve: (A) via 
(^) and Lanczos quadrature or (B) using brute 
force i.e. two further exploratory simulations at 
shifted (3 and shifted k. It is straightforward to 
estimate the relative work involved. We find that, 
to achieve the same absolute error, method (A) 
requires less than 5% of the work for (B) |p^ . 

We have conducted tests of the formalism start¬ 
ing with an ensemble of 40 configurations of a 
8^ X 24 lattice at (/3 i,ki) = (5.2,0.1340) gener¬ 
ated via Hybrid Monte Garlo. In particular, we 
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Figure 1. Predictions for the shift 5(3 from ref¬ 
erence point (5.2,0.1340) required to keep the 
specified observables constant when k is changed 
by 6k = —0.0005. Circles correspond to Wilson 
loops (1x1, 1 X 2, 2 X 4); diamonds to C 7 r(t) 
{t = 0, 1, 2); square to tq. 

Table 2 

Test of matching for the average plaquette < 
P >. The additional error in the matched value 
is systematic, reflecting the uncertainty in the 
matched beta value: 5(3 = 0.0199(10). 


(/^? Csw j 

<P> 

(5.200,2.0171,0.1340) 

(5.220,1.9936,0.1335) 

0.5286(3) 

0.5290(3)(4) 


have estimated the shift SI3f required to compen¬ 
sate for a change Sk = —0.0005 while the value 
of F is held fixed. Figure shows a compila¬ 
tion of shifts SPf for a variety of observables F: 
Wilson loops, the Sommer scale vq and pion 
correlators C^(t). In Table || we show that the 
shift prediction for the average plaquette (1x1 
Wilson loop) is accurate. An independent HMC 
simulation has been carried out at the predicted 
matched point (5.220,0.1335), using the appro¬ 
priate value of Csw There is good agreement be¬ 
tween the matched values of < P >. The partial 
derivatives (||,^) of the average plaquette mea- 
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Figure 2. Static potential on an 8^ x 24 lattice 
at (5.2,0.1340) and (5.220,0.1335) where the av¬ 
erage plaquette values match. 


sured at the reference point are 


d < P > 
Ok 


7.5(17), 


d<P> 


0.182(40) 


In Figure |^, we see some evidence that the 
shifts required to match observables which are 
sensitive to longer range physics are somewhat 
larger (perhaps 50%). With the modest num¬ 
ber of test configurations used, the statistical ev¬ 
idence is not strong. However, Figure conhrms 
this point. One sees that matching the plaquette 
ensures matching of the short distance potential 
but not the slopes of the potentials, and hence tq. 
Also shown in the hgure is the potential with a 
shift in K but not (3. There is no matching at any 
distance in that case. 


4. OTHER APPLICATIONS 

The Lanczos quadrature technique allows one 
to build relatively efficient stochastic estima¬ 
tors (j|) of the fermion determinant suitable for 
a variety of purposes including: 

• matching observables (see above) 

• tuning approximate actions ||] 


• tuning exact algorithms based on these 

• choosing parallel tempering ensembles [|lT|. 

In particular, a suitably truncated version may be 
used to model the low eigenmodes of the deter¬ 
minant fl^ while using other gauge invariant op¬ 
erators (e.g. Wilson loops) [|| to account for the 
remainder. Preliminary tests |l^ have demon¬ 
strated that relatively few Lanczos iterations are 
required to obtain a good quantitative account 
of the fermion determinant. Both the truncation 
parameters and the gauge loop parameters must 
be tuned, using the above techniques, to obtain 
a useful representation of the true QCD action 
and to ensure a good acceptance of configurations 
proposed by the gauge part of the update. Such 
approximate actions can of course be made exact 
via an accept/reject step, given sufficiently accu¬ 
rate tuning j|]. 
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